Localization in light nuclei 
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We investigate the presence of spatial localization in nuclei using a method that maps the nucleon 
same-spin pair probability and is based on the density-matrix. The method is used to study spatial 
localization of light nuclei within the Hartree-Fock approximation. We show that the method pro- 
vides an alternative tool for studying spatial localization in comparison to the localization observed 
from maxima in the nuclear mass density. 
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INTRODUCTION 



Clustering phenomena in light nuclei have always been 
an intriguing aspect of nuclear structure physics. The- 
oretical understanding of why and how conglomeration 
of nucleons to subunits within a nucleus results in an in- 
crease in stability remains an actively investigated ques- 
tion. In particular, alpha clustering in light nuclei has a 
long history [l|-[4| and suggests the existence of configura- 
tions resembling the formation of nuclear molecules [5|-l7| . 
It has also been suggested that neutron rich isotopes of 
some light nuclei may give rise to new types of cluster 
structures @, 13 ■ Most of the theoretical analyses for the 
cluster structures have been performed with the a pri- 
ori initialization in terms of clusters and effective inter- 
actions, which are determined such as to reproduce the 
binding energies and scattering phase shifts of these con- 
figurations. On the other hand, nuclear structure calcu- 
lations based on the independent-particle approximation 
or density functionals also manifest cluster-like substruc- 
tures as marked concentration of density in the visual- 
ization of the total nuclear mass density. For example, 
Hartree-Fock (HF) calculations for light nuclei often show 
such formations [9j, however since the HF single-particle 
states are generally spread across the whole nucleus they 
are delocalized, which makes the entanglement of these 
substructures in terms of the single-particle orbitals very 
difficult. Furthermore, the identification of cluster and 
shell structures based only on the mass density may be 
an oversimplification since it is missing other aspects of 
the many-body system, for example the kinetic energy 
density or density gradients, which may help to provide a 
more detailed understanding of the underlying structure. 
Finally, with the rising popularity of the density func- 
tional approach in nuclear physics it may be desirable to 
have a new localization measure that stems directly from 
the nuclear density-matrix, since all of the information is 
contained in this quantity. 



II. THE LOCALIZATION MEASURE 

A. Outline of formalism 

An alternative measure of localization had been devel- 
oped in the context of a mean-field description for elec- 
tronic systems [l(| • A fermionic mean-field state is fully 
characterized by the one-body density-matrix 

P ^(r,r') = ^0 Q (ra)C(rV). (1) 

The probability of finding two nucleons with the same 
spin at spatial locations r and r' (same-spin pair proba- 
bility) for isospin q is given by 

P qa {v,v') = Pqa {v)p qa {v>) - \ Pqaa {v,v')\ 2 , (2) 

where p gC r(r) = /9 gC rcr(r, r) is the local density. The con- 
ditional probability for finding a nucleon at r' when we 
know with certainty that another nucleon with the same 
spin and isospin is at r is 

^- ■ (3) 
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Since we are interested in the localization aspects of this 
probability it is sufficient to consider only the local short- 
range behavior of the conditional probability, which one 
can obtain by performing a spherical averaging over a 
shell of radius d about the point r and then Taylor ex- 
panding the resulting expression to get (lfjj 
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where T qa and j 9(T arc the kinetic energy density and cur- 
rent density given by 
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The reason for writing Vp g(T explicitly is to emphasize 
that to have a smooth behavior of the quantities calcu- 
lated below it is essential to calculate all quantities di- 
rectly from the wavefunctions. The expression shown in 
Eq. ((4]) suggests the definition of a localization measure 
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Light nuclei are often expected to contain pro- 
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which is also valid for time-dependent Slater determi- 
nants fill ]. It is important to remember that D qa is the 
short-range limit of the conditional like-spin pair proba- 
bility and may contain correlations that are not evident 
in simple one-body obscrvables, such as the mass density. 
The localization measure defined by Eq. JSJ) is a reverse 
relation, e.g. the larger the probability of finding two 
like-spin particles in vicinity of each other the smaller 
the value of D. For this reason it is customary to define 
a reversed and normalized localization measure 



Mr) 



TF 
'qa 



TqaPqa ~ | [V ' Pqa? j 



qa 



Pqo 



r TF 



(0) 



'qa 



where t™ is the Thomas-Fermi kinetic energy density. 



The latter is used to provide a natural scale which then 
allows to define a dimensionless measure. The current 
density vanishes in the static case which we will consider 
in the following. 



nounced a-particle sub-structures. Such a sub-structure 
means that in a certain region of space only an a par- 
ticle is found which in turn is signaled by C = 1 in this 
region. In fact, an a sub-structure is a correlation of 
four particles: p t, p 4-> n t> an d n X- Thus it is signaled 
only if we find simultaneously for all four corresponding 
localization functions C qa ~ 1. In the following, we will 
consider mainly N = Z nuclei for which the four dif- 
ferent particles have very similar wavefunctions. In this 
case, it suffices to consider, pars pro toto, only one local- 
ization function. Furthermore, it should be noted that 
a full identification of a-cluster sub-structures requires 
also to check the correlations between the four nucleons 
gathering in a "localized" region of the nucleus. The lo- 
calization function is just the first step to identify those 
regions, namely the minimum necessary condition. 



III. RESULTS AND DISCUSSION 

In our calculations, the static HF equations are solved 
on a Cartesian three-dimensional mesh without any sym- 
metry assumptions. The grid spacing was 1 fm with a box 
size of (—15.5, +15.5) fm in each dimension. The Skyrme 
energy functional was employed with the parametrization 
SkI3 [13[ . The spatial derivatives are calculated using the 
fast Fourier transform and periodic boundary conditions 
are employed, except for the Coulomb potential, which 
is calculated with boundary conditions at infinity as de- 
scribed in Rcf. flUl. 



B. Limiting cases and interpretation 

This criterion ([6]) is known in electronic systems as 
electron localization function (ELF) and it is used as one 
ingredient to analyze the bond structure of molecules in 
the static fTp| and dynamic domain [ll| . The information 
content of the localization function is understood from 
considering limiting cases. 

The extreme case of ideal metallic bonding is realized 
for homogeneous matter where r = tJJ' . This yields 
C = 5, a value which thus signals a region with a nearly 
homogeneous Fermi gas as it is typical for metal elec- 
trons, nuclear matter, or neutron stars. The opposite 
regime are space regions where exactly one single-particle 
wavefunction of type qa contributes. This is called lo- 
calization in molecular physics. Such a situation yields 
D qa (r) = 0, since it is not possible to find another like- 
spin state in the vicinity, and consequently C = 1, the 
value which signals localization. It should be noted that 
the localization function is invariant under unitary trans- 
formations am ong st the single-particle wavefunctions in 
a Slater state [12|- In the nuclear case, it is the a par- 
ticle which is perfectly localized in this sense, i.e. which 
has C = 1 everywhere for all states. Well bound nuclei 
show usually metallic bonding and predominantly have 



A. Ground states of N = Z nuclei 

Fig. Q] shows an x-z-cut of the localization function © 
for even-even N = Z nuclei from A = 4 to A = 20. The 
left panel shows the proton localization criterion C p -\ com- 
plemented in the right panel by the corresponding total 
density. As mentioned above, the states are spin sym- 
metric which yields identical localization plots for spin- 
up and spin-down. Moreover, for light N = Z nuclei 
proton and neutron localizations are very similar due to 
the small Coulomb interaction. (For neutron rich iso- 
topes this is no longer true as we will show below.) The 
color (gray scale) coding is shown on top of each column 
and remains the same throughout the column. The posi- 
tion of the density contour at half nuclear matter density 
(p = 0.08 fm~ 3 ) is indicated with color cyan in the maps 
of proton localization. One should keep in mind that the 
maxima and minima of the total nuclear density need 
not be correlated with that of the localization function, 
which is a topological quantity to describe localization 
(see also Fig. [5] and discussion thereof). The top panel 
of Fig. [T] shows the calculations for the 4 He nucleus. As 
we have described previously we see a perfect localization 
with C = 1 in all relevant regions where p > 10~ 4 fm -3 . 
Smaller densities lead to erroneous results for C due to the 
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proton localization total density p 
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FIG. 1: (Color online) Color map (gray scale) plots of proton 
localization (left column) and total density in fm -3 (right 
column), for the Z = N nuclei up to 20 Ne. The position of 
the density contour at half nuclear matter density (p = 0.08 
fm~ 3 ) is indicated with color cyan (light gray) in the maps of 
proton localization. 



very subtle cancellations required. The strongly prolate 
8 Be shows very distinct localization pattern with perfect 
localization in the left and right halves of the contour 
plane and much smaller localization in the contact re- 
gion where the wavefunctions overlap. As can be seen 
this is much more pronounced in comparison to the to- 
tal mass density plot. Here, it is probably reasonable to 
conclude that 8 Bc can be considered as an a-a molecule. 



With this version of the Skyrme force the ground state 
of 12 C is oblate deformed as shown in the right pane of 
Fig. [1] One may be tempted to consider this as a pla- 
nar arrangement of three a particles. A slight indication 
of that may be spotted in the localization plot. But it is 
not well developed, the configuration is too compact, and 
shows preferably metallic binding as we can see from the 
localization (left column) which stays safely in a regime 
C ~ 1/2. The strongly bound 16 O nucleus mostly shows 
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FIG. 2: (Color online) Density profile and localization func- 
tion for the 16 O nucleus. 



a localization value of C = | throughout as one would 
have expected. Its density is known to have a dip at the 
center [l5(. This cannot be discriminated in the density 
plot here but can be observed as a region of lower lo- 
calization in the localization map plot. To examine this 
further we have repeated the same calculation for 16 
using the SLy4 interaction [l6| . In Fig. [3] we show a cut 
through the profile of the density and similarly through 
the localization function. We observe that the central 
dip in the total density is barely visible. The localization 
function, however, shows a very pronounced dip indicat- 
ing a strong and irreducible overlap of all wavefunctions 
in this center region. Note, furthermore, that the max- 
ima of mass density and localization do not coincide. The 
localization has a preference towards the surface where 
the lower density enhances the chance of finding one pre- 
vailing wavefunction. 

Finally, the last panel of Fig. [T] shows results for the 
strongly prolate 20 Ne nucleus. The localization map 
shows two regions of high localization at the outer ends 
and a ring of somewhat enhanced localization at the cen- 
ter around the elongation axis. One can interpret this 
as a quasi-molecular a- 12 C-a configuration. The a sub- 
structures on both sides are almost as well developed as in 
8 Be. We have also computed the further series of N = Z 
nuclei, 24 Mg, 28 Si, 32 S, 36 Ar, and 40 Ca. These nuclei 
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are increasingly compact and all show basically metallic 
binding similar to 12 C and 16 shown here. 



B. Strongly deformed configurations of light N — Z 

nuclei 
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FIG. 3: (Color online) As figure [T] but for chain- like isomers 
of 12 C and 16 Q. 



larger wavefunction overlap. One interesting point about 
the 12 C linear-chain configuration localization plot is that 
in studying the dynamical formation of this chain state, 
as it was done in Ref. [l7j . we have observed that the 
dynamical vibrations of the mass density resembled the 
localization plot with only the equilibrium shape hav- 
ing the triple-a structure. This is consistent with the 
kinetic interpretation of the localization function, sug- 
gesting that kinetic energies of the same-spin pairs peak 
mostly around the ends of the linear-chain. 

neutron localization total density p 
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Very light N = Z nuclei are likely to display iso- 
meric, or at least transiently stable, configurations which 
are very elongated and resemble chains of alpha parti- 
cles [2Q|. For somewhat heavier N = Z one finds often 
shape coexistence with strongly prolate deformed nuclear 
configurations [21j . Such less bound and spatially more 
extended configurations are more likely to allow for a 
sub-structures. We thus have also considered such iso- 
meric configurations for a number of light N = Z nuclei. 
These configurations were found numerically by starting 
the static iteration from a sufficiently prolate configura- 
tion such that the iteration converged to the elongated 
isomeric state. Chain configurations were found immedi- 
ately for 12 C and 16 while the heavier systems preferred 
to maintain a compact core between the a satellites. It is 
to be noted that these configurations are stable minima 
in a mean field calculation. They may hybridize with the 
ground state in correlated calculations. Still such con- 
figurations may show up as transient configurations in 
nuclear reactions [20j . 

Fig. [3] shows the total density and localization plots 
for the linear-chain states of 12 C and 16 nuclei. For 
both the density suggests an a-chain structure which is, 
indeed, corroborated by the localization that also shows 
three or four clearly separated maxima, C ~ 1. The 
region of high localization is very large at both ends, 
but much smaller for the maxima in the interior due to 



FIG. 4: (Color online) As figure [T] but for stretched isomers 
of 24 Mg and 28 Si. 

Fig. |4] shows strongly prolate (not yet chain-like iso- 
mers which lie higher in energy) isomers of 24 Mg and 28 Si. 
Unlike the compact ground-state configurations these iso- 
mers indicate interesting molecular substructures. One 
may interpret 24 Mg as a a- 12 C-a-a molecule and 28 Si as 
a-a- 12 C-a-a. Again, the outermost a's are best devel- 
oped with large regions of high localization. The inner 
a's have already degraded localization due to neighboring 
wavefunctions from both sides. 

C. An example for N > Z: The 20 C chain 

Recently, much interest has been devoted to the study 
of cluster co nfig urations for neutron-rich isotopes of light 
nuclei JfJ [H, Il8| . In particular the linear-chain configu- 
rations of C isotopes and their stability against bending 
modes has been of interest. For nuclei with N > Z where 
proton and neutron wavefunctions are naturally differ- 
ent the search for a sub-structure requires a simultane- 
ous analysis of proton and neutron localization. To that 
end we consider also as a localization the combination 
\/CpfC n f . The spin-up and spin-down wavefunctions are 
still degenerate such that it suffices to consider one of the 
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FIG. 5: Color map plots of localizations and density for the 
linear-chain configuration of 20 C. Lower left: proton localiza- 
tion. Upper left: neutron localization. Lower right: a local- 
ization (= y/C p ^C n ^). Upper right: total density (in fm -3 ). 
The position of the density contour at half nuclear matter 
density (p = 0.08 fm -3 ) is indicated with color cyan (light 
gray) in the maps of localization. 



spins. In Fig. Owe show proton, neutron, and a localiza- 
tion plots for the linear-chain isomer of the 20 C nucleus. 
As expected, due to the neutron excess of 20 C the local- 
ization plots for neutrons and protons look considerably 
different. The protons show more distinct regions with 
high localization value in comparison to the neutron case, 
where the wavefunctions have more overlap due to the 
large number of neutrons. The a localization is the obvi- 
ous average of the two left panels. In spite of the neutron 
cloud from the excess neutrons, there appears still some 



faint a sub-structure at the edges of the chain. It is also 
interesting to observe that the total mass density does 
not show any pronounced features due to the smoothing 
effect of the surplus neutrons while the localization plots 
still reveal noteworthy structures. 



IV. CONCLUSION 

In summary, we have applied a localization measure 
which was developed originally for analyzing bonding 
structures in molecules to a study of alpha sub-structures 
in light nuclei. The localization function is obtained di- 
rectly from the density-matrix, in the mean-field approxi- 
mation. It depends on kinetic-energy density and current 
density, in addition to the mass density. It can be easily 
implemented for density functional theory calculations of 
nuclear structure. One of the fundamental reasons why 
the new localization measure is such an excellent predic- 
tor of correlation and localization is due to the fact that 
it incorporates the kinetic energy of the relative motion 
of spin-parallel nucleons at a particular point in space 
in addition to the mass density for the system In 
most cases this localization function shows more detailed 
localization or clustering features in comparison to the 
total mass density. Results for N = Z nuclei up to 40 Ca 
show that pronounced localization, associated with a- 
particle substructures, appear only for the strongly pro- 
late ground states of 8 Be, 20 Ne, and of course trivially for 
4 He. All other nuclei are more compact and show metal- 
lic binding. However, stretched isomers of light nuclei 
often show convincing a structures, particularly well de- 
veloped for the a chains of 12 C and 16 O, but also for the 
prolate 24 Mg and 28 Si isomers. In the future we also plan 
to study the new localization function in time-dependent 
HF calculations of systems related to nuclear molecular 
configurations. 
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